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Abstract 



Baym, Blaizot, and Zinn-Justin have recently used the large approxima- 
tion to calculate the effect of interactions on the transition temperature of 
dilute Bose gases. We extend their calculation to next-to-leading-order in 
1/N and find a relatively small correction of -26% to the leading-order result. 
This suggests that the large approximation works surprisingly well in this 
application. 

I. INTRODUCTION 

Second-order phase transitions have universal behavior, associated with long wavelength 
fluctuations, for which critical exponents and other universal quantities can often be success- 
fully calculated using renormalization group techniques. For most such systems, the short 
distance physics is hopelessly complicated. In contrast, the phase transition of a dilute, 
interacting Bose gas provides a fascinating example where physics becomes simpler, and 
perturbative, at (relatively) small distance scales. For this system, it should be possible to 
marry techniques for treating long-distance critical fluctuations to a perturbative treatment 
of short distance physics, and so compute non-universal characteristics of the phase transi- 
tion. A simple example of such a non-universal quantity is the phase transition temperature 
Tc, and the effect of interactions on Tc has been explored by several authors f^-^], with a 
wide variety of theoretical results. In particular, the transition temperature has recently 
been calculated by Baym, Blaizot, and Zinn-Justin |^ in the large N approximation. For 
simplicity, they implicitly focus on the case of Bose gases with a single spin state, where 
the low-energy cross-section for atomic collisions can be parametrized by a single scattering 
length, a. As will be briefly reviewed below, the problem is flrst reduced to a calculation in 
a three dimensional 0(2) scalar field theory at its critical point. Replacing that by an O(A^) 
theory with N = 2, they find 
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in the dilute limit, where n is the number density,^] and Tq is the transition temperature of 
a non-interacting Bose gas, 
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Their result of AT^/To = (Tc — Tq)/Tq ~ 2.33 an^/^ is in good agreement with recent^ 
numerical simulations that give ATc/To ^ (2.2 ± 0.2) an^/^ The result IS surprismg 
because it seems to work much better than the large expansion of critical exponents. For 
example, for 0{N) theory, the susceptibility critical exponent 7 is P]0 
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This is not a marvelous expansion for N = 2, for which the actual value is 7 ^ 1.32. 
In this paper, we calculate the 0{N^^) relative correction to ( p..l|) . We find 
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Setting N = 2, this is only an 26% correction to the leading large N result for ATc/Tq. 
We now have ATc/Tq ~ 1.71 an^^^. Though this does not agree as well with the quoted 



1 For simplicity, we consider a uniform Bose gas, where n is fixed. Alternatively, in an arbitrarily wide 
harmonic trap, n should be interpreted as the actual density at the center of the trap at the transition 
temperature. 

^ As alluded to earlier, there have been several different theoretical results and simulation results obtained 
by various methods (e.g. ||l|-^), giving a large range of values for the coefficient of an^^^ in ATc/Tq. There 
has also been some experimental data on the ^He-Vycor system |^] , which superficially seems to fit well an 

early theoretical estimate of Stoof §], which is ATc/Tq ~ (167r/3) C{3/2)-'^/^an^/^ ~ 4.66 an^/^. However, 
the detailed interpretation of this data is unclear. In that experiment, the Helium atoms are confined to an 
interconnected network of channels in the porous Vycor glass, and, for the low-density data that appears 
to fit Stoof, the interpaticle spacing is the same order of magnitude as the widths of the channels. Ref. [|j 
simply assumes that the system can be modeled by a free Bose gas with (i) an effective mass for the atoms 
that is extracted experimentally, but (ii) the same scattering length as for bulk Helium, which is moreover 
taken from theoretical modeling. Because of these assumptions, the apparent agreement with Ref. j^] should 
be treated with caution. 



For Bose gases, a more physical example of a critical exponent is = f — 0.540(2/A^) — 0.470(2/A^)^ + 
0{N^^), whose actual value is ~ 0.67 for N = 2. The fact that 0(2) critical exponents should be identified 
with Bose gas exponents is not completely trivial. A uniform, non-relativistic Bose gas is a constrained 

system: the particle density n is fixed. This constraint causes the critical exponents x = (a,/3,7, t') of 
the actual system to be related |f^ to the standard exponents x = (a,/?, 7,:^) of the field theory by (i) 
a = —a/{l — a), and x = x/{l — a) for the others, if a > 0, or (ii) 5; = a; if a < 0. The actual value of 
a for the 0(2) model is believed to be — 0.007 ± 0.006 |ll[ . If negative, there is no difference between the 
exponents; if positive, there is in principle a very tiny olrference. This relation explains, by the way, the 
difference between mean- field theory exponents for the 0(2) model (e.g. a — 1/2) and the exponents of a 
non-interacting Bose gas (e.g. a = —1). 
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simulation result, the moderatley small size of the correction supports the proposition that 
the large N expansion works surprisingly well for T^. 

In the remainder of this introduction, we review the long-distance 0(2) effective theory 
for Bose condensation and then review the arguments of about how to calculate ATc/T^. 
In Sec. |T|, we review the leading-order calculation in large N as done in [0. In Sec. |ITT|, we 
go on to calculate the next order in An appendix explains how to calculate some of 

the basic 3-dimensional integrals that appear at that order. 



A. Review of effective theory 

The basic assumption throughout will be that the average separation n~^^'^ of atoms is 
large compared to the scattering length a. This can also be expressed as X{Tq) ^ a, where 
A is the thermal wavelength 

It is well known that, at distance scales large compared to the scattering length a, an appro- 
priate effective theory for a dilute Bose gas is the second-quantized Schrodinger equation, 
together with a chemical potential fi that couples to particle number density and a 
contact interaction that reproduces low-energy scattering. The corresponding Lagrangian is 

£ = ^* (-ih dt-—V''-fx]ij + iri^Y- (1-6) 

y 2m J m 

In this context the corresponding mean-field equation of motion is called the Gross-Pit aevskii 
equation.^ As with any effective theory, there are corrections represented by higher- 
dimensional, irrelevant interactions (in the sense of the renormalization group), such as 
(t/'*?/^)^ and ilj*'V^ip. However, higher and higher dimension operators are parametrically less 
and less important if the distance scales of interest are large compared to the characteristic 
scales (a) of the atomic interactions. The term in the Lagrangian (|1.6|) is in fact the 

lowest-dimension irrelevant interaction, and it is adequate for computing the leading-order 
effects of interactions in the diluteness expansion.^ 

Now treat the system at finite temperature using Euclidean time formalism. The field 
can then be decomposed into frequency modes with Matsubara frequencies Un = 2Txnk^T /%. 
At sufficiently large distance scales (^ A), and small chemical potential (|yu| <^ k-^T), the 
— (/;,^/2m)V^ — terms in ( |1.6| ) become small compared to the 0{hujn) time derivative 



For a review, see |p^ . 

^ At short distances, the dt and terms of the action j dtdPx C determine that times scales as (length)^ 
and the scahng dimension of ijj is (length) 

^ For a discussion of analyzing corrections in this language, see ref. jlSl , which extended earlier work on 
corrections by refs. ImI. A similar discussion for Fermi gases may be found in ref. psf. 
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term, provided ri 7^ 0. The non-zero Matsubara frequency modes then decouple from the 
dynamics, leaving behind an effective theory of only the zero-frequency modes ipQ. Roughly, 
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with f3 = l/ksT. In detail, the parameters of (|1.7] ) are renormalized by coupling to the non- 
zero modes, and there are again corrections in the form of irrelevant (and even marginal) 
interactions. However, these effects are all suppressed in the dilute limitQ and do not affect 
the computation of AT^/Tq at leading order in an^^^. 

It is then convenient to write ipo = h~ {mkBTY^'^{(j)i + i(j)2) so that the effective action 
S = H/T becomes a conventionally normalized 0(2) field theory: 
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where (j) is understood to be a 2-component real vector (( 
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B. Review of ATJTq 

Our effective theory depends on two as yet undetermined parameters — r and u, or equiv- 
alently /i and T. One constraint comes from fixing particle number density n: 

n={r^) = ^{<P'). (1.10a) 

At the critical temperature, a second requirement is that the system have infinite correlation 
length, which requires 

=r + n(0) =0, (1.10b) 

where 11 (p) is the proper self-energy of the field. The two equations (|1.10|) determine the 
two unknowns r and u, and hence Tc. As noted by Baym et ai, the density equation ( |1.10a] ) 
can be rewritten as 



^ For the 3-dimensional effective theory (1.7), the short-distance scahng dimension of -00 is (length) 
the (V^qi/jq)^ interaction is relevant, and a (i/'o'^/'o)'^ interaction would be marginal. Even though marginal, 
this last interaction can be ignored at the order of interest in the diluteness expansion because it has a 
small coefhcient. For example, conside r the term that would arise directly from the prese nce o f a correction 
g3{ip*tp)^ to the original Lagrangian (1.6). That would lead to a g3(V'oV'o)'^ term in (1.7) which, after 
rescaling, would become a term proportional to {mgs/h \'^){(fP)'^ in (1.8). Since A ~ ^{Tq) oc rt^^/'^ at the 
transition, the coefhcient of this term is high order in the diluteness expansion in n^^^. Similarly, an effective 
(ipQipo)^ term arising from the 4-point interactions {ip*ip)'^ and from integrating out physics at the scale A 
(due, for example, to non-zero Matsubara modes) would give rise to a ((^^)^ term in (1.8) with coefficient 
proportional to u^\^ cx A~'^. 
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where the last equahty uses ( |1.10b|) . Throughout this paper, we will use the notational 
shorthand 
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for momentum integrals. (Technically, p is a wave number rather than a momentum, but 
we will use conventional fi=l nomenclature, even though we have not set to 1.) 

The expression (|1 . 1 1|) for the density is ultraviolet (UV) divergent and so receives con- 
tributions from short distance scales where the effective theory breaks down. This could be 
handled by appropriately regulating the effective theory and then perturbatively correcting 
the UV contribution. As pointed out by Baym et ai, it is simpler to instead consider the 
difference n — no(T), where no(T) is the same expression in the absence of interactions (i.e., 
with n set to zero): 
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no{T) represents the density a non- interacting Bose gas has if its transition temperature is 

am 



T. It is given by inverting (1.2): 
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This formula cannot be derived directly in the effective theory (|1.8|) , but the difference n — rio 
in ( |1.13| ) is insensitive to the UV and so can be. 

The above constraints are entirely adequate to systematically determine AT^ in the large 
N expansion, but there is a convenient way to simplify the bookkeeping a bit. Baym et al. 
give a simple argument that, to leading order in the density expansion. 
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where the factor of 2/3 in ( |1.13| ) arises from the relation T oc n^^^ . Combining (|1.15| ) with 



we can summarize as 
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to leading order in an^^^ . It's also useful to rephrase this, again in terms of the fields of 
the effective theory, as 
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I 

FIG. 1. Diagrams contributing to A{(fP) at leading order in 
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FIG. 2. Bubble chains. Unbroken lines denote flavor index contractions. 

where 

A(02) = (02) - (02)n_o. (1.18) 



Note that the problem of calculating A(02) from the action ( |1.8| ), subject to the con- 
straint (|1.10b|) , has only one dimensionful scale in it: u. The length scale of this problem, 



which will be the length scale of the physics that determines AT^/Tq, is therefore 

\2 

(1.19) 

a 

by dimensional analysis. In the dilute limit A(Tc) ^ a, this length scale is large compared 
to A, which justifies use of the 0(2) effective theory ( p. .81) . 



II. REVIEW OF LEADING ORDER IN 1/N 



We now review the leading large- calculation of (0^) — (02)n^o; and hence of ATc/Tq, 
by Baym et al. The details of our calculation are slightly different than theirs, and we will 
introduce techniques needed to proceed to higher order. We start with the standard large- iV 
generalization of the 0(2) scalar field theory ( |1.8| ) to an 0{N) scalar theory: replace by 
an A^-component vector and treat Nu as fixed in the N (yo limit. The reader should keep 
in mind that u is therefore order 1/N. Standard N power counting of Feynman diagrams 
consists of a power of m ~ 1/N for each 4-point vertex and a power of N for each flavor 
trace. 

The set of diagrams that determine A{(jP) at leading order in 1/A^ is depicted in Fig. 
ID, where the dashed line denotes bubble chains, as shown in Fig. H. (For comparison, the 
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FIG. 3. Diagram representing the non-interacting result (0^)n^o- 



diagram for (0^)n^o is shown in Fig. |^.) The cross denotes an insertion of the operator 
0^, whose expectation we are computing. There is a simple way to summarize the effect 



on diagrammatic perturbation theory of the rcj)^ term in the action ( |1.8| ) and the constraint 
( iLlObI) that r = -n(0). 

Rule 1: Use massless (gapless) scalar propagators l/p^ when evaluating diagrams, 
ignoring the term in the action. But whenever there is a one-particle irre- 
ducible sub-diagram X that represents a contribution to the proper self-energy 
Il{p), ther^ replace X{p) by X{p) — X(0). 

We note for later reference that, for the purpose of this rule, a diagram that is cut in two 
pieces only by cutting a single internal dashed line is still one-particle irreducible, because 
cutting the bubble chain represented by a dashed line corresponds to cutting two (f) lines 
(Fig. I). 

The bubble chain sum shown in fig. is given by 



= -N~^Fp, (2.1) 
where 

Fp = WfTT ^2.2a) 
and So(p) represents the basic massless bubble integral (not summed over flavors) 

In d = 3 dimensions. 



^ This rule is unambiguous for calculating expectations such as {(j)^)- It is potentially ambiguous for 
calculating the free energy — for example, a diagram like Fig. ^ but without the cross on it. In that case, 
it is ambiguous which sub-diagrams would be considered self-energy insertions. A systematic way to treat 
the perturbation theory in all cases is to treat the r(f? term in the action as a perturbation, include it in 
Feynman diagrams as a 2-point vertex, and then set r = —11(0) order by order in perturbation theory. 
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Putting everything together, the diagram of fig. IT] gives 



A(<^^) = - / 
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As pointed out by Baym et al. in [|^, the above integral is not absolutely convergent 
in three dimensions, and one must be careful to consistently regulate the theory before 
proceeding. Integrals that are not absolutely convergent are at best ambiguous — they depend 
on the order one chooses to do the integrations. For example, if one evaluates ( |2.4| ) directly 
in three dimensions, doing the angular integrations first, then the I integration, and then the 
p integration, the result is zero. This is not in fact the correct answer. We will discuss this 
issue in some detail in order to justify the correctness of our procedure for later evaluating 
higher-order diagrams. 

Baym et al.^s preferred method for the leading-order calculation is to use dimensional 
regularization and evaluate everything in c? = 3 — e dimensions. This is difficult at next order 
in the loop integrals we shall encounter are sufficiently complicated that evaluation in 
d = 3 — e dimensions seems hard. Our strategy will be to instead always reduce diagrams to 
well-defined three-dimensional integrals, which are simpler to evaluate. We imagine starting 
with some consistent regularization scheme, like dimensional regularization, and will now 
discuss how to manipulate the integrals so that they will be absolutely convergent if we set 
d = 3. We assume in what follows that the UV regulator respects parity and is invariant 
under shifts p ^ p + k of loop momenta p. 



Let's look at the divergences that cause absolute convergence of the integral (2^) to fail 
in three dimensions. A simple one to correct is the behavior for 1 fixed and p —>■ oo. The large 
p piece of the p integration then behaves as /p p ■ which is logarithmically UV divergent 
(from the point of view of absolute convergence). This can be remedied by rewriting the 
regulated version of (|2.4|) by using the freedom to change the integration variable p to — p: 



LO 



ip /4 
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Pi 



211 



Pi 
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Now, if we throw away the UV regulator, the 1 fixed, large p divergence is gone. This 
sort of divergence is trivial, easy to remedy, and won't have much practical impact on our 
calculations (given the order in which we will eventually do integrations). We'll simply 
acknowledge the issue in later calculations, without emphasizing it, by writing ( p.5|) as 



LO 



ip 



II 



Pi 
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where the subscript ± means that one should average the expression with p — > — p (or 
equivalently with 1 — —1). 

Unfortunately, even (|2.5|) is not absolutely convergent. There is still a logarithmic UV 
divergence associated with I and p simultaneously becoming large {I ^ p ^ oo), as can be 
seen by simple power counting and the fact that Fp approaches a non-zero constant for large 
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p. Return to considering (|2] 
divergence by rewriting 



with a UV regulator still in place. We can eliminate the UV 
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The second integral vanishes, as can be seen by changing integration variable p 
its first term. So 
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A(0')lo 
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This is now UV convergent because Fp — ^ as p ^ oo. However, we have traded 
the logarithmic UV divergence for a logarithmic infrared (IR) divergence, associated with 



V 
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We now need some sort of infrared regulator. One physically motivated possibility for 
consistently regulating the infrared would be to consider the system infinitesimally above 
the critical temperature, so that all the massless scalar propagators should be replaced 
by massive ones l/(p^ + M^), where the mass M represents the inverse correlation length 
This defines an absolutely convergent integral in 3 dimensions, and the limit M — 
would be taken only after the integrations. 

Massless propagators will be much easier to deal with, however, in higher-order 
calculations. As a practical matter for computing diagrams, we prefer to introduce as few 
massive propagators as possible. It would be convenient, for example, to IR regulate 
by introducing M only in the 1/P propagators: 
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M-»0 



Ip {p + M^y 



|i + p| 



(2.9) 



J ± 



where Fp is still defined in terms of the massless bubble integral, as in (|2.2| ). One might worry 
that an ad hoc procedure of putting masses only on some propagators could be inconsistent, 
so let us argue more carefully. Return to the UV regulated version of (|2.5|) and note that the 
integral is not sensitive to the the region of integration where / is infinitesimal, because this 
particular integral is IR convergent. There's then no reason we can't modify the infrared 
behavior of the integrand for infinitesimal I, without affecting the integral. So, for instance. 
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But now, again rewriting Fp = {Fp — F^o) + F^o, the same steps as before reproduce (|27 

Now that we have an absolutely convergent integral ( p. 91) , we can do the integration in 
three dimensions and in any order we choose. It's convenient to do the 1 integral first: 



1 (/2 + M2)2 
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The "±" prescription makes no difference to this particular integral, because the 1 integration 
by itself is completely convergent without it. Note that naively setting M to zero at this 
stage would give the incorrect, zero result mentioned earlier. Instead, we have 
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The overall factor of M in the numerator is canceled by a linear IR divergence in the p 
integration, which is cut off by M. 

For small M, the integral ( p.l2|) is dominated^ by p ~ M. So, in the limit of M — 0, we 
can simplify the calculation slightly by replacing — F^q by Fq — F^. So 

A(0^)lo = (Fo - F^) g^^2(pf+M2) = = • ^^-^^^ 

When combined with the formula ( p..l7| ) for ATc, this reproduce Baym et a/.'s leading large 
N result ( p. .11) , in which N has been set to 2. 



III. NEXT ORDER IN 1/iV 



The diagrams which contribute to A(0^) at next order in are shown in Figs. ^ and 
H. The diagrammatic expansion comes from the standard introduction of an auxiliary field 
a, represented by the dashed lines.0 The 0{N) action of (|1.8|) is rewritten as 
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(fx 



-|A(/.|2 + ir02 + i0V-— (T^ 
2 ' ^' 2 ^ 2 ^ 6u 



(3.1) 



The a propagator is then turned into the bubble chain of fig. |^ by resumming the basic 
massless bubble of fig. ^ into the a propagator. Technically, this is accomplished by trivially 
rewriting (p.l|) as 



S — Sq -\- iSsubtractions + J ^ 0^(7 , 



2 C^-pl 
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subtractions 



k r(l)-p(f)p + \ cx_p A^So(p) CTp 



(3.2a) 
(3.2b) 
(3.2c) 



with Fp and So given by (p.2|). The terms designated 5'subtractions may be ignored if one 
follows the previous Rule 1 as well as 



^ Some readers raay worry that the integral ( ^.12 ) is dominated by arbitrarily small p ^ M — ^ 0. They may 
worry because at sufficiently small momentum our perturbative propagators are no longer good approxima- 
tions to the full propagators. The full scalar propagators, for example, actually scale like l/l'^^'^ rather than 
l/l"^ at small I (<C Nu), where the critical exponent 77 is 0{N~^). The differen ce be comes significant when 
I ^ Nuexp{—ri) ~ iVuexp[— 0(l/iV)]. One might worry that the sensitivity of (|2.12| ) to p ^ is a sign that 
naive large N perturbation theory must break down. It is important to realize, mthe present case, that this 
infrared sensitivity is simply an artifact of our mathematical manipulations on the infrared-safe expression 
( ^.9|) . Regardless of whether one used some sort of infrared- improved propagator in (2^), that expression 
is not sensitive to far-infrared momenta. It is sensitive to momenta Nu, for which there is nothing wrong 
with a large- TV expansion based on perturbative propagators. 
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For a very quick review of standard large N, see, f or e xample, section 2.1 of chapter 8 of 



Some 

people might prefer to replace a by ia in the action (3J^), so that the Euclidean path integraTfor a is 
convergent, but it matters not at all for the purpose of large N perturbation theory. 
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(a) 



(b) 



(c) 




(d) 




(e) 



text. 



FIG. 4. Next-to-leading order diagrams for A{(jP). p,q, 1, 1' label loop momenta, as used in the main 




+ 




FIG. 5. Diagrams for the a self-energy T,{p) at 0{N''^). 



Rule 2: Do not include any diagrams that have the one-loop bubble, fig. as a 
sub-diagram. 

Note that Rule 1 eliminates any tadpole sub- diagrams, such as fig. |^. Formal large 
counting of diagrams is simply to count a factor of for each a propagator and a factor 
of for each (p loop. The important momentum scale of the problem will be the scale 
p ~ Nu = 0{N'^), where the a propagator (|2.1|) makes the transition from its small p 
behavior {Fp oc p) to its large p behavior {Fp constant). Some authors like to completely 
integrate out of ( |3.1| ), but we prefer to retain it, as there is then a more transparent 
relationship between Feynman diagrams and the corresponding Feynman integrals. 

In evaluating the diagrams of fig. |^, we shall borrow techniques from ref. [|17|, where 
somewhat related diagrams were evaluated in gauge theories with large numbers of scalars. 



FIG. 6. The one-loop bubble diagram. 



11 



FIG. 7. Example of a tadpole diagram. 



Our strategy will be to do the loop integrals first, and then tackle the remaining integrals 
associated with a propagators. 



A. Diagram a 



Let's start with Fig. Ha. The corresponding integral is 



pq 



|l + Pl 



|l + q| 



(3.3) 



J ± 



As written, this integral is absolutely convergent and can be evaluated, without regulariza- 
tion, directly in three dimensions. To do the 1 integration, however, we find it convenient to 
temporarily introduce an IR regulator mass M. We may then separately integrate each of 
the terms of the integrand, which are not individually IR convergent. We can also use M as 
a trick for reducing powers of /~2_ Specifically, we rewrite the 1 integral as the M ^ limit 
of 



1 



1 {P + M2)3 
1 <P 



1 



|l + Pl 



1 

pi 



1 



1 



2 rf(M2)2 7i (Z2 + M2) 

1 (p 



Jl + q| 

1 

|T+iF 



1 

p2 



J ± 



|l + q| 



J ± 



2 d{M 



2^2 



Jl (p, q; M) - p-'Ji (q; M) - q-'Ji (p; M) + p^'q 



where 



/i(p,q;M) = 
Ji(p;M)^ 



1 



(/2 + M2)|l + p|2|I + q|2^ 

1 



1 (/2 + M2)|I + p 



2' 



au' 



(3.4) 

(3.5) 
(3.6) 



The integral + Af^)^^ is —M/4tt plus an Af -independent UV divergence, and derivatives of 
the latter vanish. It's of course not necessary to introduce J^{P + M^)~^ and this spurious UV divergence; 
one could simply evaluate J-^{P + M"^)^^ direct ly. But we find it convenient to consolidate the treatment of 
such integrals with that of the other terms in (3.4). 
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1 1 /• 1 I (P f M\ 



- 



(3.7) 



The integral Ji is straightforward to evaluate. A particularly simple way to evaluate Ji 
is to make a conformal transformation which reduces it to the form of Ji. The results of 
both integrals, and the conformal transformation between them, are discussed in Appendix 
^ All we need here are the small M expansions of those results, which turn out to be 

/i(p, q; M) 



8pg|p — q| 47rp2g2 8p^g^|p — q| 



^ M3(/ + 4p-q + g^) ^ M^(3(p ■ q)^ - pV) ^ 

127rp^g^ 16p^q^\p — q| ' 



1 



Ji(p;M) = - 



M 



M3 



+ 0(M^). 



(3.9) 



Putting everything together. 



l(P + ]\ 



{P + M2)3 

p ■ q 



|l + p|2 -2 



1 

p2 



+ 



3(p ■ a 



± 

2-1 



l + q|2 g2 



± 



SirMp'^q'^ 16p^g'^|p — q| 
+0,M,. 

lop'^g'^lp — q| 



+ 0(M) 



J ± 



J ± 



We can now set M = 0. All that will matter in the integral (|3.3| ) is the average (■ 
the angle between p and q, which is 



1 

1/6 



|l + p|2 p2 



J ± 



1 



|l + q|' 



1 



AOp^p^ 



where 

We are left with 



p> = max(p, q), p< = min(p, q). 



1 



F F 



N Jpd 40p» p< 27r^N Jo 



oo rp F F 

p^dp / q^dq ^ 







(3.10) 
■)e over 

(3.11) 

(3.12) 
(3.13) 



The remaining integrals are easy to do, with the result 



u 



7r2 5 



157r4 \ 6 4 



(3.14) 
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B. Diagram b 



Fig. corresponds to 



-1 ' F F B 
pq 



(3.15a) 



5pq = / - 



i/M |i + p| 



|l + p + q|2 g2 



p4 



1 



1 



|p + q|2 g2 



(3.15b) 



This contribution to A (02) is again absolutely convergent if the subscript ± is taken to mean 
averaging over p — > — p and also over q — > — q. It is convenient to now rewrite the I integral 
as the Ml, limit of 



1 



1 (/2 + M2)2 |(|l + p|2 + M2)2 

d d 



1 



1 



|l + p + q|2 g2 



(p2 + M^y 



1 



|p + q|2 



dm) dm) - ^^2(p;M„M,) 



(p2 + Ml) 



|p + q|2 g2 



Ml 
An 



where 



J2(p, q-M,,M2) = 



/2(|I + p|2 + Mi2)(|l + q|2 + Ml 



2\ ' 



(3.16) 



(3.17) 



J2(p;Mi,M2) 



l(/2 + Mi2)(|l + p|2 + M|)- 



(3.18) 



The results for I2 and J2, and their small Mi, M2 expansions, are given in Appendix 0. The 
final result for the I integration, after taking the Mi, M2 — > limit, is 



B 



pq 



q - 2p(p ■ q) - 3g(p • q)^ 
8p2g2|p + q|3 



(3.19) 



with angular average 



(^pq)6 



d{p - q) 

Ap^q 



(3.20) 



where 6{p — q) is the step function (1 for p > q; for p < q). The remaining integrals are 
easy to do, giving 



U / TT 



37rM 6 A ' 



(3.21) 
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C. Diagram c 



Fig. ^ can be evaluated as the others, but the final integrals are a bit more complex. 
The diagram gives 



(3.22a) 



1 



|1 + p|^|l + qp|I + p + qp p2g2|p + qp 
The 1 integration can be performed using methods similar to before: 



J ± 



a 



d 



pq 



lim 



iJ(p,q;M) 



1 



p^g^lp + qp 



M 
47r 



(3.22b) 



(3.23) 



ff(p,q; M)^ 



1 



(/2 + M2)|l + p|2|l + q|2|l + p + q|2 



(3.24) 



H can be reduced to the basic integrals Ji and Ji encountered previously by rewriting the 
numerator 1 in ( p.24| ) as 



1 



(Z2 + M^) + |1 + p + q|2 - |1 + p|2 - |1 + q|' 



2p ■ q + M2 

and then expanding the integrand into the corresponding four terms: 
1 



(3.25) 



ff(p,q;M) 



2p ■ q + M2 



[/i(p, q; 0) + /i(p, q; M) - /i(p + q, q; M) - /i(p + q, p; M)] . 



Using the expansion ( p.8|) of Ji, one obtains 



pq 



1 + (p ■ a 



,-2 



|p-q| |p + q| 



p q - ip qj 

8p2 g2 |p _|_ q| 



J ± 



This simphfies, after applying the [. . .]± prescription, to 



p ■ q - IP ■ qj 



pq 



16 p2 g2 



|p + q| |p - qr 



Angular averaging yields 

(Gpq) 

where 



8pl,a;2(l + a;2) 



X + 



Sinh a; 



(3.26) 



(3.27) 



(3.28) 



(3.29) 
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X = P</p>- 



(3.30) 



The remaining integrals over p and q are no longer so trivial. Notice first that for an 
arbitrary function f{x) one can rewrite 

8m /"i , x^lnx 



F F 



N Jpd p> 
so that 



ri p p 



u 



1 , xlnx 
ax- 



X 



+ 



37r4 Jo (1 - x) 
Sinh~^x 



37r4 7o (1 - x) 
This is easiest to evaluate numerically, giving 

where c ~ 0.463715. We also have an analytic result:^ 



l + x2 (l + a;2)3/2 



cu 



TT 



48 



1 + ln(l + V2) 



3 ^(3,X8). 



where 



[S^Xs)- gs gs + ys + g. 13^ + 15^+17^ 



is a particular case of Dirichlet's L-function, and 

L(3,X8) = 0.958380454563 



/(x), (3.31) 

(3.32) 

(3.33) 

(3.34) 

(3.35a) 
(3.35b) 



D. Diagram d 

Fig. ^d corresponds to 



D 



pq 



1 



1' r|l' + p|2|l' + q|2 7i /4 



|l + p|2|l + q|2 p2g2 



/i(p,q, M)-p-'q''{- 



±1 

M 

47r 



p • q 

g4p3g3 1 p — q 



2 ■ 



(3.36a) 



(3.36b) 



Our inelegant, brute force method for obtaining this result is borrowed from a footnote of ref. |rp]. The 

hard part is the Sinh~ term. We change variables from x io y — x -\- yl + aP . This turns Sinh^ x into 
Iny. \nx can be written as a sum of terms of the form ln(?; — a), and the change of integration variable 
makes the rest of the integrand a rational function of y. We then split this rational function apart by partial 
fractions and do each int egr al, yielding di- and tri-logarithms of various arguments. Finally, we use a zoo of 
polylogarithm identities |18| to simplify the answer. 
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Here the subscript ±1 means we implicitly average over 1 
Doing the remaining integrals by brute force, we find 



u 



7_ 
12 



C(3) 



TT 



6 



-1 for absolute convergence. 



(3.37) 



where ( is the Riemann zeta function. It's interesting to note that 7r^/6 can also be written 
as C(2). 



E. Diagram e 

The final class of diagrams, fig. correspond to the leading order diagram, fig. with 
the replacement 



(3.38) 



where S(p) represents the contribution to the a self-energy at next-to- leading order, shown 
in the diagrams of fig. ||. Making this substitution in the leading-order calculation (|2.13|) 
gives 



^0 ~ ■^a 

327r2 



(3.39) 



The nice feature of this relation is that we only need to calculate S (p) in the small and large 
p limits. 

The first self-energy diagram in fig. |^ contributes 



Si(p) 



1 



i /-. ^"^^P+ql /2|i^ + p|2|i^ + p + q|2 4 i2\i^ + p|2|i2 + p + q|2 • 



1 



(3.40) 



For small p, the integration is dominated by /i, ^2 ~ P and q ~ Nu. So we can ignore li, I2, 
and p compared to q and write|^ 



1 2 



Nu 



2^2 ■ 



(3.41) 



This diagram has a quadratic IR divergence for p = 0. In contrast, the other two diagrams 
only diverge as linear x log when p = 0, and so behave as p~^lnp for small p. In summary. 



S(P) 



Nu 
487r2p2 



+ 0(p"Mnp) 



(3.42) 



for small p. One may check by power counting diagrams and sub-diagrams that S(oo) = 0. 
Then 



13 For general p, the result is Nu T.i{p) = 2471^^ z-^Rc[2 Li2(-z) - 2 Li2(2 + z) + \tt% where z = 48p/Nu, 
and Li2{z) = — dx ln(l — x)/x is the dilogarithm function. One may double check that the p — > limit 
agrees with ( [3.4l| ). 
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^0 



16 

3^ 



ail' 



d0 



T 



U 

"6^ 



(3.43) 



(3.44) 



We conclude by mentioning one technical subtlety, glossed over above, concerning abso- 
lute convergence. The integration corresponding to the substitution ( |3.38| ) in the leading 
order analysis is 



T 



|i + p| 



(3.45) 



In contrast to the analogous leading-order expression ( ^.6] ), this integral is not absolutely 
convergent in the infrared (/ ~ p ^ 0), though it is convergent in the UV. One might 
therefore worry about the ad hoc introduction of an IR regulator M in the calculation of 
this graph. However, this worry is easily bypassed by rewriting 



ip 



(jFp — JFq) 



II 



Pl 



(3.46) 



± 



which should be understood as regulated in the UV. The UV-regulated integral of the jFg 
factor vanishes. Eq. ( |3.46|) is now convergent in the IR, but logarithmically divergent in the 
UV, just as the original leading-order integral ( |2.6| ) was. One can now follow through the 
same argument as in the leading-order case to introduce an IR regulator and then remove 
the UV divergence, where Fp in the leading-order analysis in now replaced by J-'p — JFq. The 
result is still ( p.39|) . 



Summary 



Summing all the diagrams then yields the total NLO contribution: 



NLO 



U 



I ^(3) -I- iLvr^ + ^ ln(l + V2)-- L(3, Xg) 
4^^ ^ 2 240 48^2 ^ 3 ^ ' 



with L{?),x%) given by (|3.35| ). Combining with the leading-order result ( 2.13 ), 



A(0^)nlo _ 0.527198 
which is the relative NLO correction we presented for in ( |L4| 



(3.47) 



(3.48) 
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One may also check this answer by direct, brute- force calculation of all the diagrams associated with fig. 
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APPENDIX A: BASIC INTEGRALS 

Let's begin with the integral J2{p', Mi, M2). This is quite easy to do by standard methods 
(for example, by introducing a Feynman parameter), and gives 

J2(p;Mi,M2) 



1 (/2 + M2)(|l + p|2 + M|) 



Svrp \{Mi + M2y+p' 

Sp 47rp2 ^ I27rp4 ^ V ; V ; 

The integral Ji(p; m) of ( p.6|) is simply the special case Ji(p; M) = J2(p; M, 0). 

The integral /2(p, q; Mi, M2) of ( p.l7| ) can be related to J2(p, q; Mi, M2) by generalizing 
a trick presented in ref. '^\- The idea is to change integration variables from I to its 
conformal inversion 1 = 1//^. The integration measure changes as (27r)^^(i^Z = [2Ti)~^l~^d^l. 
Propagators can be written in terms of the new variable I as 

I = (A2) 



where 



K^-T^- (A4) 



Making this change of variables. 
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/2(p,q;Mi,M2) 



(p2 + M2)(g2 + M|) 



1 + 



P 



X 



p2 + Mf 



+ 



Ml 



g2 + M| 



+ 



Mo 



1 



(p2 + M2)(g2 + M|) 



J2 



P 



q 



+ M|^ 
Ml 



Mo 



• (A5) 



For the application of this paper, the relevant terms in the small Mi, M2 expansion of I2 are 

1 1 



/2(p,q;Mi,M2) 



8pg|p — q| 47r|p — qp y 



^ + ^]+0{Ml) + 0{Ml) 



M1M2 /Ml M2 



47r|p — q|^ \ (f 



+ 



+ 0(Mf) + 0{Ml) 



+MlMl 



[3pg - 2(p^ + g^)(p ■ g) +pg(p ■ g)^ 

8p2g2|p — q|5 



0(M5). (A6) 



The integral Ji of (|3.5|) is related by /i(p, g; M) = hip, p — g; M, 0) and gives 

Cos"i(2a;2^-l) 



/i(p,q;M) 



where 



87rM|p-g|2Jcu-2-l 



M|p - g| 



pq 



'(p2 + M2)(g2 + M2) 
The small M expansion is given in (|3.8|) . 



(A7) 



(AS) 
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